perm filename MIX.F4[MUS,LCS] blob
sn#058113 filedate 1974-01-08 generic text, type T, neo UTF8
00005 C LOAD WITH MIXIO AND UNPACK
00010 IMPLICIT INTEGER (A-B,X-Z)
00100 DIMENSION A(1024),B(1024),X(3),Y(3),Z(3)
00110 EQUIVALENCE(X1,X(1)),(X2,X(2)),(X3,X(3))
00120 1,(Y1,Y(1)),(Y2,Y(2)),(Y3,Y(3))
00130 1,(Z1,Z(1)),(Z2,Z(2)),(Z3,Z(3))
00150 TYPE 10
00160 ACCEPT 11,Q,ITOTL
00170 IF(Q.EQ.0)Q=2000
00200 N='MUSAA'
00300 CALL GETFIL(N)
00400 N='INPUT'
00500 CALL GETFI2(N)
00510 4 N='XXXXX'
00555 CALL PUTFI2(N)
00600 N=1
00700 1 CALL FASTIN(A(1),1024)
00800 CALL FASTI2(B(1),1024)
01000 DO 2 K=1,1024
01100 CALL UNPACK(A(K),X1,X2,X3)
01200 CALL UNPACK(B(K),Y1,Y2,Y3)
01300 DO 3 M=1,3
01310 G=Y(M)
01320 IF(G.GT.2047)G=G-4096
01330 CC60 G=(G+2048)/Q
01335 60 G=G/Q
01340 6 H=X(M)
01350 IF(H.GT.2047.)H=H-4096.
01400 3 Z(M)=G*H
01450 C MULITPLIES INPUT*MUSAA IN 20:70 RATIO.
01500 2 A(K)=Z1*16777216+Z2*4096+Z3
01510 CALL FASTO2(A(1),1024)
01600 5 IF(N.GE.ITOTL)GO TO 40
01700 N=N+1024
01800 GO TO 1
02200 40 CALL FINFI2
02210 10 FORMAT(' TYPE INPUT MAXAMP',$)
02220 11 FORMAT(F,I)
02300 END